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ABSTRACT 



Aims. The purpose of the present study is to research the origin of planetary companions by using a precise radial velocity (RV) 
survey. 

Methods. The high-resolution spectroscopy of the fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at Bohyunsan 
Optical Astronomy Observatory (BOAO) is used from September 2008 to June 2012. 

Results. We report the detection of two exoplanets in orbit around HD 208527 and HD 220074 exhibiting periodic variations in RV 
of 875.5 ± 5.8 and 672.1 ± 3.7 days. The RV variations are not apparently related to the surface inhomogeneities and a Keplerian 
motion of the planetary companion is the most likely explanation. Assuming possible stellar masses of 1.6 + 0.4 and 1.2 ± 0.3 M Q , we 
obtain the minimum masses for the exoplanets of 9.9 +1.7 and 11.1 ± 1.8 Mj up around HD 208527 and HD 220074 with an orbital 
semi-major axis of 2.1 ± 0.2 and 1.6 ± 0.1 AU and an eccentricity of 0.08 and 0.14, respectively. We also find that the previously 
known spectral classification of HD 208527 and HD 220074 was in error: Our new estimation of stellar parameters suggest that both 
HD 208527 and HD 220074 are M giants. Therefore, HD 208527 and HD 220074 are so far the first candidate M giants to harbor a 
planetary companion. 

Key words, stars: planetary systems - stars: individual: HD 208527: HD 220074 - stars: giant - technique: radial velocity 



1. Introduction 

Most exoplanets have so far been detected by using spectro- 
scopic radial velocity (RV) and photometric transit methods at 
the ratio of two to one. While the transit method is limited by 
using strongly biased towards short-period companions with a 
period as short as a few days, the RV technique is relatively 
less restricted, it can detect planetary systems with a period ap- 
proximately from a day to ~ 38 years (i.e. 47 Uma described in 
Gregory & Fischer 2010). 

It has been demonstrated that the RV method has a limited 
range of stellar spectral types of stars chosen as targets. Planet 
searches generally focus on late F, G, and K dwarfs because 
they are bright enough to obtain a high signal-to-noise ratio 
(S/N) in the high-resolution spectra and have an ample num- 
ber of spectral lines for the RV measurements. Furthermore, 
finding planets around the Sun-like stars is a big issue in ev- 
ery aspect. Unfortunately, there are no more candidate stars in 
the main-sequence (MS) stage for ground-based observations. 
Recognizing the importance of exoplanet formation and its evo- 
lution, several exoplanet survey groups have begun exploring 
evolved G, K giants (Frink et al. 2002; Setiawan et al. 2003; Sato 
et al. 2003; Hatzes et al. 2005; Dolinger et al. 2007; Han et al. 
2010) and low-mass (M dwarf) stars (Endl et al. 2003; Bonfils et 
al. 2011). 

A giant star has sharper spectral lines than a dwarf star, 
which allow more precise measurement of RV shift. RV surveys 
of giant stars are also interesting because we can detect planets 



* Based on observations made with the BOES instrument on the 1.8- 
m telescope at Bohyunsan Optical Astronomy Observatory in Korea. 



of stars that are considerably more massive than the Sun. For 
MS stars of the same mass this would be difficult because they 
rotate fast. Compared to MS stars, however, giants show blended 
RV variations due to various surface variations in the line pro- 
files, such as chromospheric activities, and stellar pulsations. It 
thus makes it difficult for the giant star to purify the RV origin 
through the examination process. Also, it is not so easy to deter- 
mine the masses of the host stars, because stars of different mass 
are located in roughly the same region of the H-R Diagram. And, 
the interpretation of the data is not easy, because only massive 
planets with relatively long orbital periods have been detected 
for giant stars, whereas surveys of MS stars often find low-mass 
planets with short orbital periods. 

Of the exoplanets detected by the RV method, just over 12% 
are discovered around giants. However, no planet has been dis- 
covered around M giant stars because they are rare or maybe 
because planet formation is weak around M giants unlike M 
dwarfs, which may have a significant probability (up to 13%) 
of harboring short-period giant planets (Endl et al. 2003). Short- 
term oscillations on M giants are rarely observed, yet they show 
periods in excess of one day (Koen & Laney 2000). Three- 
dimensional simulations (Freytag & Hofner 2008) predict typ- 
ical photometric periods of about several tens to a few hundred 
days, similar to those of RV variations in K giants. 

For the past four years, we have conducted precise RV mea- 
surements of 40 K dwarfs. Here, we present two stars with long- 
period and low-amplitude RV variations. In Sect. 2, we describe 
the observations and data reduction. In Sect. 3, the stellar char- 
acteristics of the host stars are reinterpreted, including chromo- 
spheric activities and Hipparcos photometric variations. The RV 
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Table 1. Stellar parameters for the stars analyzed in the present paper. 



Parameters 


HD 208527 


HD 220074 


Reference 


Spectral type 


K5 V 


Kl V 


Hipparcos (ESA 1997) 




- 


M2 III 


Kidger et al. (2003); Tsvetkov et al. (2008) 




Ml III 


M2 III 


Derived" 


m„ [mag] 


6.4 


6.4 


Hipparcos (ESA 1997) 


M v [mag] 


- 1.24 


- 1.52 


Famaey et al. (2005) 


B-V[mag] 


1.68 


1.66 


Hipparcos (ESA 1997) 




1.70 + 0.02 


1.68 ± 0.01 


van Leeuwen (2007) 


V-K [mag] 


4.14 


4.39 


Cutri et al. (2003) 


age [Gyr] 


2.0+ 1.3 


4.5 ± 2.8 


Derived* 


Distance [pc] 


320.2 ±61.1 


290.2 ± 45.5 


Famaey et al. (2005) 


RV [km s- 1 ] 


4.79 ± 0.06 


-36.89 ±0.21 


Famaey et al. (2005) 


Parallax [mas] 


2.48 ± 0.38 


3.08 ± 0.43 


van Leeuwen (2007) 


T tS [K] 


3950 


5080 


Wright et al. (2003) 




4643 ™ 


4971 +741 

y ' -643 


Ammons et al. (2006) 




4035 ± 65 


3935 ± liO 


This work 


[Fe/H] 


-0.09 + 0.16 


- 0.25 ± 0.25 


This work 


logg 


1.6 + 0.3 


1.3 ±0.5 


This work 




1.4 + 0.2 


1.1 ±0.2 


Derived* 


R* [Re] 


51.1 ±8.3 


49.7 ± 9.5 


Derived* 


Mi, [M Q ] 


1.6 + 0.4 


1.2 ±0.3 


Derived* 


v ro , sin [km s -1 ] 


3.6 


3.0 


This work 


P rot / sin i [days] 


660.0 ± 116.6 


832.4 ± 156.7 


Derived" 


Vmicro [kms -1 ] 


1.7 ±0.3 


1.6 ±0.3 


This work 



<a) See text. <b) Derived using the online tool (http://stevoapd.inaf.it/cgi-bin/param). 

measurements and orbital solutions will be presented in Sect. 4. 
Finally, in Sect. 5, we discuss our results from our study. 

2. Observations and reduction 

Since 2003, we have conducted a precise RV survey for ~ 300 
stars, including 55 K giants (Han et al. 2010), ~ 200 G giants 
(Omiya et al. 2008), 10 M giants, and 40 K dwarfs. The obser- 
vations were carried out as part of the study on the 40 K dwarfs. 
We used a fiber-fed high-resolution (R = 45 000) Bohyunsan 
Observatory Echelle Spectrograph (BOES; Kim et al. 2007). 
The BOES spectra covered a wavelength region from 3500 to 
10 500 A. In order to provide precise RV measurements of the 
K dwarfs, an iodine absorption (I2) cell was used with a wave- 
length region of 4900 - 6000 A. Each estimated S/N ratio at an 
I2 region was about 200 with a typical exposure time ranging 
from 10 to 20 minutes. The long-term stability of the BOES was 
demonstrated by observing the RV standard star t Ceti, which 
was constant with an rms scatter of 6.8 m s _1 over the time span 
of our observations (Lee et al. 201 1, 2012). 

We obtained 24 and 32 spectra each for HD 208527 and HD 
220074, respectively, from September 2008 to June 2012. The 
basic reduction of spectra was performed using the IRAF (Tody 
1986) software package and DECH (Galazutdinov 1992) code. 
The precise RV measurements related to the I2 analyses were un- 
dertaken using the RVI2CELL (Han et al. 2007), which is based 
onamethodby Butler et al. (1996) andValenti et al. (1995). The 
RV measurements for HD 208527 and HD 220074 are listed in 
Tables|2]and [3]and shown in Figure [4] and [7] 

3. Stellar characteristics 

3. 1 . Fundamental parameters 

The atmospheric parameters were determined based on the 
TGVIT (Takeda et al. 2005) code. To estimate the parameters, 




3955 3960 3965 3970 3975 3980 

Wavelength (A) 



Fig. 1. The Ca II H line cores for samples. The vertical dotted 
line indicates the center of the Ca II H regions. There are slight 
central emissions that exist in the center of the Ca II H core fea- 
ture for HD 208527 and HD 220074, which are ambiguous to 
confirm the samples to be active. 61 Cyg A is shown for com- 
parison. 



we used 169 (HD 208527) and 145 (HD 220074) equivalent 
width (EW) measurements of Fe I and Fe II lines. They resulted 
in r eff = 4035 + 65 K, log g = 1.6 + 0.3, v micro = 1.7 + 0.3 km s" 1 
and [Fe/H] = - 0.09 + 0.16 for HD 208527, and r eff = 3935 ± 
110K, logg= 1.3 ±0.5, Vmicro = 1.6 + 0.3 km s- 1 and [Fe/H] = 
- 0.25 + 0.25 for HD 220074. The projected rotational velocities 
were derived, v ro t sin i = 3.6 and 3.0 km s , from the widths of 
the spectral lines regarding the line-broadening functions. Based 
on the rotational velocities and the stellar radius of 51.1 + 8.3 
and 49.7 + 9.5 R , we derived the range of the upper limit of 
the rotational period of 660.0 + 1 16.6 (HD 208527) and 832.4 + 
156.7 (HD 220074) days. 
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8658 8660 8662 8664 8666 

Wavelength (A) 



Fig. 2. Line variations in the Ca II 8662 A spectral region. While the active star 61 Cyg A shows the variations of EW, there are 
invisible variations for HD 208527 and HD 220074 at the Ca II 8662 A core feature. They show typical line profiles of an inactive 
star. 



Stellar masses were estimated from the theoretical stellar 
isochrones by using their position in the color-magnitude dia- 
gram based on the ones by Bertelli et al. (1994) and Girardi et 
al. (2000). We also adopted a version of the Bayesian estimation 
method (J0rgensen & Lindegren 2005; da Silva et al. 2006) by 
using the determined values for T e s, [Fe/H], M v , and parallax. 
With our estimated parameters, we yielded a stellar mass of 1.6 
+ 0.4 (HD 208527) and 1.2 + 0.3 M (HD 220074). The basic 
stellar parameters are summarized in Table [1] 

HD 208527 ( = HR 8372 = HIP 108296) has been known to 
be a MS star with a spectral type of K5 V (ESA 1997) and an 
apparent magnitude of 6.4. The absolute magnitude of - 1.24, 
however, indicates that the star could be located at a giant stage 
in the H-R diagram. The values of the surface gravity and stellar 
radius also show that it may be a giant star. The surface gravity 
of the star is consistent with a giant but inconsistent with a dwarf. 
Stars that have the same spectral class but a different luminosity 
class have quite different optical versus infrared colours. Bessel 
& Brett (1988) provide the intrinsic V-K colours for dwarfs and 
giants. The V-K colour of 4.14 is close to an Ml giant (V-K = 
4.05) rather than a K5 giant (V-K = 3.60) according to Bessel 
& Brett (1988). Thus, HD 208527 should be classified as an M 
giant rather than as a K dwarf star. 

The classification of HD 220074 ( = HR 8881 = HIP 
115218) has caused controversy. Whereas the result of the 
Hipparcos photometric observation estimated that the star re- 
mains in a MS stage with a spectral type of Kl V, recent reports 
showed a spectral type of M2 III (Kidger et al. 2003; Tsvetkov 
et al. 2008). The colour index of HD 220074 was estimated to be 
approximately 1 .7 by Hipparcos data, which was reduced afresh 
by Leeuwen (2007) and shows the same result. It means that the 
star should be a spectral type of a late K or early M. However, 
the effective temperatures have led to confusion over the spectral 
classification of HD 220074. Studies have shown it to be ~ 5000 
K (Wright et al. 2003; Ammons et al. 2006), which indicates 
the spectral type of an early K and it agrees with the Hipparcos 



classification of a Kl-type. On the contrary, our estimation of 
an effective temperature (T s h = 3935 K) shows that the spectral 
type is close to a late K or an early M. HD 220074 also could 
be classified as a giant star according to the value of the absolute 
magnitude M v = - 1.52, as well as a low surface gravity and a 
large stellar radius compared to that of the dwarf star. The V-K 
colour of 4.39 is similar to an M2 giant (V-K = 4.30) in Bessel & 
Brett (1988), which also indicates that the star is likely to be a gi- 
ant. As a result, HD 220074 should be an M2 giant star (Kidger 
et al. 2003; Tsvetkov et al. 2008). 

3.2. Chromospheric activities 

The EW variations of the Ca II H & K lines are frequently used 
as chromospheric activity indicators because they are sensitive 
to stellar activity and the variations could have a decisive effect 
on the RV variations. The emissions in the Ca II H & K core are 
formed in the chromosphere and show a typical central reversal 
in the presence of chromospheric activity (Pasquini et al. 1988; 
Saar & Donahue 1997). FigureQ~]shows individual features in the 
Ca II H line region. Due to a low S/N of the blue-most region, 
the spectrum is not clear enough to resolve the emission feature 
in the Ca II H line core. It may be the result of some invisible 
emission features or scatters in the both samples. Even if there 
really are slight central emissions compared to 61 Cyg A (K5 V) 
showing atmospheric activity, which is one of the samples of 40 
K dwarfs, the features are too weak to confirm the existence of 
core emission features categorically. 

Chromospheric activity also has been measured by using the 
Ca II triplet lines (Larson et al. 1993; Hatzes et al. 2003; Han et 
al. 2010). Among the triplet lines (8498, 8542, and 8662 A) , Ca 
II 8662 A is uncontaminated by atmospheric lines near the core. 
Changes in the core flux of the Ca II 8662 A reflect variations 
in stellar chromospheric activity and qualitatively related to the 
variations in the Ca II H & K flux (Larson et al. 1993). The line 
profiles in the Ca II 8662 A spectral region are shown in Figure|2] 
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Table 2. RV measurements for HD 208527 from September 
2008 to June 20 1 2 using the BOES . 



600 900 
Period (days) 



Fig. 3. The Hipparcos photometric variations for HD 208527 
(top panel) and HD 220074 (bottom panel), respectively. The 
horizontal line indicates a FAP threshold of 1%. 



The small boxes are magnified regions of 1 A centered on the Ca 
II 8662 A to show the EW variations near the core. The compar- 
ison of 6 1 Cyg A, known as the active star, was obtained during 
the same season at the epoch (JD-2454726.198, 2454726.220, 
2454752.124, 2454833.894, and 2455018.172). The spectra of 
61 Cyg A clearly show variations of the central core, which is 
very similar to the result in Larson et al. (1993). On the other 
hand, the line features for HD 208527 and HD 220074 display 
no such variations. 

3.3. Photometric variations 

We analyzed the Hipparcos photometry data to search for pos- 
sible brightness variations, which may be caused by the rota- 
tional modulation of cool stellar spots. The available photom- 
etry database originates from 88 Hipparcos measurements for 
HD 208527 (from November 1989 to December 1992) and 121 
Hipparcos measurements for HD 220074 (from January 1990 
to December 1992). They maintained a photometric stability 
down to rms scatters of 0.010 (HD 208527) and 0.012 (HD 
220074) magnitude, which correspond to 0. 15% and 0. 18% vari- 
ations over the time span of the observations. Figure [3] shows 
the Lomb-Scargle periodograms (Lomb 1976; Scargle 1982) of 
the Hipparcos data for both stars. A rather large peak appears 
at about 301 days with a false alarm probability (FAP) of 0.2 + 
0.1% (HD 208527). A FAP is a metric to express the significance 
of a period by a bootstrap randomization technique (Press et al. 
1992; Kiirster et al. 1999). But, there are no additional signifi- 
cant peaks at other periods in both periodograms. 



4. Radial velocity variations and origin 

4.1. HD 208527 

Period analysis was applied to appraise the significance of pe- 
riodic trends or variabilities of the time series. A Lomb-Scargle 
periodogram, which is a useful tool to investigate long-period 
variations for unequally spaced data, was used to calculate the 
RV time series. The time series of our RV measurements for 
HD 208527 are listed in Table [2] and shown in Figure |4] (middle 
panel). Figure [5] shows a significant power for HD 208527 at a 
period of 875.5 days (top panel) with a FAP of less than 10~ 5 and 
the residual after subtracting the superior period (bottom panel). 



JD 


ARV 




JD 


ARV 




-2 450 000 


m s 


m s 


-2 450 000 


ms~' 


ms~' 


4736.119458 


-91.0 


16.5 


5454.126840 


87.2 


9.3 


4752.152971 


-84.2 


11.3 


5697.263528 


-148.0 


20.2 


4833.916842 


-176.6 


8.0 


5698.288703 


-98.7 


10.8 


5018.194241 


-82.3 


11.0 


5711.267459 


-125.7 


12.8 


5083.979529 


-3.6 


8.6 


5820.129124 


-80.6 


9.4 


5131.919794 


-5.7 


10.9 


5842.035628 


-208.3 


11.1 


5169.956708 


69.6 


11.1 


5911.006597 


-138.1 


10.7 


5183.885600 


53.1 


10.2 


6024.330791 


51.7 


12.1 


5224.889144 


129.5 


19.0 


6026.309597 


38.4 


15.3 


5311.277212 


79.7 


10.4 


6069.265534 


134.2 


11.5 


5321.264398 


120.4 


11.1 


6086.197407 


202.6 


30.0 


5356.307074 


183.2 


9.2 


6088.318959 


94.0 


15.1 



100 - 
SJ 
-100 - 



300 
200 
100 


-100 
-200 
-300 

100 


-100 



a = 6.8 (m/s) 



HD 208527 




*tr= 119.5 (m/s) 
tj(O-C) = 39.3 (m/s) 



5100 5400 5700 

JD-2450000 



Fig. 4. Variations of the RV standard t Ceti (top panel), RV curve 
(middle panel), and rms scatter of the residual (bottom panel) for 
HD 208527 from September 2008 to June 2012. The solid line is 
the orbital solution with a period of 875 days and an eccentricity 
of 0.08. 



There is a highly significant power at a period of 875.5 days and 
no significant peak exists in the residual. 

HD 208527 is reclassified as a giant star in this work, rather 
than a dwarf. Thus, HD 208527 may have surface inhomo- 
geneities which are not uncommon to giant stars. Stellar rota- 
tional modulations of surface features can create variable asym- 
metries in the spectral line profiles (Queloz et al. 2001), and in 
particular, the variations in the shapes of the spectral line may 
cause the RV variations. The different RV measurements be- 
tween the high and low flux levels of the line profile are de- 
fined as a bisector velocity span (BVS). In order to calculate the 
BVS, we selected three strong lines with a high flux level and un- 
blended spectral features, which are Ni 1 6643.6, Ti 1 6743. 1, and 
Ni 1 6767.8 A. The lines are located beyond the I2 absorption re- 
gion and show a high flux level throughout the whole wavelength 
region. We estimated each BVS of the profile between two dif- 
ferent flux levels with central depth levels of 0.8 and 0.4 as the 
span points, that avoided the spectral core and wing where the 
difference of the bisector measurements are large. The BVS vari- 
ations of HD 208527 as a function of RV are shown in Figure|6] 
The slopes are 0.02 (Ni I 6643.6), 0.01 (Ti I 6743.1), and 0.03 
(Ni I 6767.8 A), respectively; There are no obvious correlations 
between the RV and BVS measurements. 
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Table 3. RV measurements for HD 220074 from September 
2008 to June 20 1 2 using the BOES . 



300 600 900 1200 1500 

Period (days) 

Fig. 5. A Lomb-Scargle periodogram of the RV measurements 
for HD 208527. The periodogram shows a significant power at a 
period of 875.5 days (top panel) and after subtracting the main 
frequency variations (bottom panel). The horizontal lines indi- 
cate FAP thresholds of 0.1%, 1%, and 10% (top to bottom). 
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19.4 


5018.233079 


-285.0 
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5911.065978 
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RV (m/s) 

Fig. 6. The examinations of the origin of the B VS variations for 
HD 208527. It shows the relations between the RV and BVS 
variations for three spectral lines (Ni I 6643.6, Ti I 6743.1, and 
Ni I 6767.8 A). The solid, dashed, and dotted lines mark the 
slopes of the three spectral lines, respectively. 
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JD-2450000 

Fig. 7. Variations of the RV standard t Ceti (top panel), RV curve 
(middle panel), and rms scatter of the residual (bottom panel) for 
HD 220074 from September 2008 to June 2012. The solid line is 
the orbital solution with a period of 672 days and an eccentricity 
of0.14. 



HD 208527 reveals a periodic signal in RV variations, and 
the best fit Keplerian orbit is found with a P = 875.5 + 5.8 days, 
a K = 155.4 + 3.2 m s"\ and an e = 0.08 ± 0.04. Figure 
shows the RV curve as a function of time for HD 208527 and the 
residual after extracting the main frequency. Assuming a stellar 
mass of 1 .6 + 0.4 M Q , we derived a minimum mass of a planetary 
companion m sin i — 9.9 + 1.7 Mj up at a distance a — 2. 1 AU from 
HD 208527. 

4.2. HD 220074 

The time series of our RV measurements for HD 220074 are 
listed in Table [3] and shown in Figure [7] A Lomb-Scargle pe- 
riodogram of the data exhibits a dominant peak at a period of 
672.1 days in Figure[8](top panel) with a FAP of less than 10~ 5 . 
Figure [8] (bottom panel) shows a Lomb-Scargle periodogram of 
the residual after subtracting the superior period and no signifi- 
cant peak exists. 

The BVS was measured by using the same spectral features 
and method as that for HD 208527. The BVS variations of HD 



220074 as a function of RV are shown in Figure|9] The slopes are 
0.10 (Ni I 6643.6), - 0.07 (Ti 1 6743.1), and - 0.01 (Ni 1 6767.8 
A), respectively. The correlation between the BVS and RV is not 
visible. However, we want to point out that this can also be due 
to too small BVS variations, as expected from Saar & Donahue 
(1997)orHatzes (2002). 

The RV variations of HD 220074 can be well-reproduced by 
Keplerian orbit with a P = 672.1 ± 3.7 days, a K = 230.8 ± 5.0 
m s , and an e — 0.14 + 0.05. The RV variations and residuals 
for HD 220074 are shown in Figure [7] Assuming a stellar mass 
of 1.2 + 0.3 M Q , we derived a minimum mass of a planetary 
companion m sin i = 11.1 ± 1.8 Mj up at a distance a — 1.6 AU 
from HD 220074. All the orbital elements are listed in Table [4] 

5. Discussion 

In order to search for and study the origin of planetary compan- 
ions, the precise RV measurements around 40 K dwarfs were ob- 
tained with BOES for four years. We find that HD 208527 and 
HD 220074 show evidences of low-amplitude and long-periodic 



5 



B.-C. Lee et al.: Planetary companions orbiting M giants HD208527 and HD 220074 



15 " HD 220074 




300 600 900 1200 1500 

Period (days) 



Fig. 8. A Lomb-Scargle periodogram of the RV measurements 
for HD 220074. The periodogram shows a significant power at a 
period of 672.1 days (top panel) and after subtracting the main 
frequency variations (bottom panel). The horizontal lines indi- 
cate FAP thresholds of 0.1%, 1%, and 10% (top to bottom), re- 
spectively. 
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Fig. 9. The examinations of the origin of the B VS variations for 
HD 220074. It shows the relations between the RV and BVS 
variations for three spectral lines (Ni I 6643.6, Ti I 6743.1, and 
Ni I 6767.8 A). The solid, dashed, and dotted lines mark the 
slopes of the three spectral lines, respectively. 

RV variations. The spectral and luminosity classes for the sam- 
ples are reclassified with the result of previous studies and our 
new estimation of the stellar parameters. HD 208527, known as 
a K dwarf, is changed to an M giant and HD 220074, known as 
a K dwarf or an M giant, is confirmed to be an M giant. Thus, 
the low-amplitude RV variations may have been caused by the 
surface phenomena as well as by orbiting planetary companions. 

We found RV periods of 875 and 672 days for HD 208527 
and HD 220074, respectively. Long-term RV variations, such as 
in hundreds of days, in evolved stars may have been caused by 
three kinds of phenomena: stellar pulsations, rotational modu- 
lations by inhomogeneous surface features, or planetary com- 
panions. Of these, RV variations due to pulsations are unlikely 
because the period of the fundamental radial mode pulsation is 
several days and the mechanism to produce very long and low- 
amplitude pulsation modes is impractical. 

The lack of activity in the stars argue in favor of the planet 
hypothesis. The Ca II H & K line and Ca II 8662 A were 



Table 4. Orbital parameters for HD 208527 b and HD 220074 b. 



Parameter 


HD 208527 


HD 220074 


Period [days] 


875.5 ± 5.8 


672.1 ± 3.7 


^periastron [JD] 


2450745.3 ± 45.7 


2451158.2 ± 31.7 


K [m s" 1 ] 


155.4 ± 3.2 


230.8 ± 5.0 


e 


0.08 ± 0.04 


0.14 ±0.05 


oj [deg] 


278.0 ± 19.5 


323.0 ± 8.5 


f(m) [M Q ] 


(3.375) x 10- 7 


(8.314) x 10~ 7 


a sin i [AU] 


(1.247) x 10~ 2 


(1.412) x 10~ 2 


o-(O-C) [m s" 1 ] 


39.3 


57.4 


m sin i [M Jup ] 


9.9 ± 1.7 


11.1 ± 1.8 


a [AU] 


2.1 +0.2 


1.6 ±0.1 



used to monitor chromospheric activities, which moreover show 
strong correlations from the results of early K stars (Thatcher 
& Robinson 1993). Because of the low S/N in the Ca II H & 
K region of the BOES spectra, the variations in the Ca II 8662 
A line profile were examined instead: they do not show any no- 
table variations. The Hipparcos photometric variations and the 
BVS also showed no correlations with the RV measurements. 
Even though there exists a moderate power at a frequency of 
301 days, as compared to 875 days seen in RV variations, in 
Hipparcos photometric data for HD 208527, this does not go 
against the general trend. All these null results from the analy- 
ses of the chromospheric activities or surface inhomogeneities 
argues against the rotational modulation of surface activities as 
the origin of observed RV variations. 

However, the existence of planetary companions can explain 
the observed RV variations. Assuming the stellar masses of 1.6 
+ 0.4 and 1.2 ± 0.3 M for HD 208527 and HD 220074, the 
minimum masses of planetary companions were estimated to be 
9.9 + 1.7 and 11.1 + 1.8 Mj up with an orbital semi-major axis of 
2.1 and 1.6 AU, respectively. If confirmed, HD 208527 and HD 
220074 become the first candidates harboring planetary compan- 
ions around M giant stars. 

After removing the RV signal of the planets, the remain- 
ing RV variations are 39.3 m s _1 (HD 208527) and 57.4 m s _1 
(HD 220074), respectively. The variations are larger than the ex- 
pected for stellar activity and the errors of the measurements. 
Additional frequencies can be considered and duly weighed be- 
cause the values are significantly larger than the RV precision 
for the RV standard star t Ceti (6.8 m s~'), as well as larger than 
the typical internal error of individual RV accuracies of ~ 13 m 
s _1 (HD 208527) and ~ 22 m s" 1 (HD 220074). A periodogram 
of the RV residual variations, however, does not show any addi- 
tional periodic signal as shown in Figure [5] (bottom panel) and 
Figure [8] (bottom panel). This means that the residual RV vari- 
ations do not indicate the presence of another planet. Hekker et 
al. (2006) have shown that the rms of RV residuals of K giants 
has a median value of 20 m s _1 and this tends to increase toward 
later spectral types. Hatzes et al. (2005) also have shown that the 
high rms of residuals (~ 51 m s~') are found in the Kl giant HD 
13189 and this shows a significant short-term RV variability on 
time scales of days that is most likely due to stellar oscillations. 
That behavior is typical for K giant stars, whether periodic or not 
(Setiawan et al. 2003; Hatzes et al. 2005; Dolinger et al. 2007; 
Lee et al. 2008; de Medeiros et al. 2009; Han et al. 2010). 
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